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Abstract 

We show that the perturbation series for quarkonia energies diverges at large orders. 
This results in a perturbative ambiguity in the energy that scales as e~^/"^ where a is 
the Bohr radius of quarkonium and A is the QCD scale parameter. This ambiguity is 
associated with a nonperturbative contribution to the energy from distances of order 
A^^ and greater. This contribution is separate from that of the gluon condensate. 



The states of quarkonia which have a size much less than the characteristic scale of 
nonperturbative QCD, A~^, can be accurately described by QCD perturbation theory, just 
as positronium is accurately described by QED perturbation theory The quarkonium 
wave function is nonzero at long distances, of order and greater, so the states of quarkonia 
are also influenced by nonperturbative physics. In this article we show that perturbation 
theory is aware of this nonperturbative effect via its large-order behavior. 

Nonperturbative effects in quarkonia have previously been paramaterized in terms of vac- 
uum condensates. At short distances, the leading such effect is due to the gluon condensate 
0, §]. We argue that this effect is separate from the nonperturbative effect we investigate, 
which is due to the nonzero wave function at long distances. 

That perturbation theory is aware of (some) nonperturbative effects via its large-order 
behaviour is familiar from the study of infrared renormalons 0, ^. Renormalons correspond 
to factorial growth of the coefficients of the perturbation series, and are associated with 
nonperturbative contributions of order {A/QY, where Q is a characteristic momentum and 
p is a positive integer. Our analysis is in the spirit of infrared renormalons. We, however,find 
perturbative coefficients that grow less rapidly than factorially, and which are associated with 
a nonperturbative contribution of order e~*^/^, rather than a power of (A/Q). 

At short distances, quarkonia can be described by the QCD analogue of the Coulomb 
potential [|l|, 

4MiZr) ,1) 
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with a coupling that depends on the distance between the quark and the antiquark. At 
distances of order A~^ and greater the coupling is large, and Eq. (|l|) is invalid. To study 
the effects of the running coupling perturbatively, we expand the evolution equation for the 
coupling as an infinite sum of terms which depend on the coupling at an arbitrary fixed 
scale n, as = Q!<j(/i). 
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is the one-loop beta function for Nf fiavors of light quarks. Inserting Eq. into Eq. (|l|) 
gives 
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with 
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Diagrammatically, the potential in Eq. (^) corresponds to the sum of vacuum-polarization 
insertions in the tree-level potential, as shown in Fig. 1. In this sense, our analysis is similar 
to that of infrared renormalons. 
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Figure 1: Quarkonium formation via the exchange of a gluon. The Coulomb potential is 
modified by the insertion of vacuum-polarization subdiagrams in the gluon propagator. 



The first term in Eq. (1) is the QCD Coulomb potential with a fixed coupling. We take 
this as the unperturbed potential. Its solution involves the usual Coulomb energies and 
wave functions. We regard the remaining terms as perturbing potentials of order {(3as)'', 
and evaluate their contribution using ordinary nonrelativistic time-independent perturbation 
theory. 

To evaluate the contribution of the k^^ perturbing potential to the energy of the state with 
principal quantum number n, orbital quantum number /, and azimuthal quantum number 
mJ3 one evaluates the matrix element 

AE^l^ =< ^nlm\V^\^nlm > ■ (6) 

For simplicity we momentarily restrict our attention to the ground state. For quarks of mass 
TTiQ, the ground-state wavefunction and energy are 

^Mr) = -^e-^/^ (7) 
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where a is the Bohr radius, 
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(The calculation for arbitrary n, I, m is neither more difficult nor more enlightening; we 
present the result later.) Inserting Eq. (|^) into Eq. (H), dividing by E, and integrating over 
angles yields 

{I3a,f drr In'' fir e-^'/" (10) 
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Making the change of variables z = fir and setting c = fia/2 gives 
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^ Since the potential is radially symmetric, the states with different azimuthal quantum numbers are 
degenerate. 
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where 
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For large fc, the integral of Eq. (|I2]) is dominated by the region where the argument of the 
exponential is a maximum. We denote this maximum by 2:05 defined by the implicit equation 

zq In 2:0 = kc . (13) 

We approximate the integral by the saddle-point method, which yields 

^ V2^c zl'^ In'^ zo e-'''/" . (14) 

Thus Eq. (|ll|) becomes 
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The series for the energy correction, whose fc*'^ term is given by Eq. (|15|), is divergent. 
For large /c, the size of the fc*'^ term is governed by the factor 

= {^Cis\^z^f e-'^l" . (16) 

The term where the series begins to diverge is obtained by finding the minimum of this 
function with respect to k (keeping in mind that zq implicitly depends on k via Eq. ([T3|) ): 
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Using Eq. (13), we find that the last two terms cancel. Solving for zo gives 



A ^ 
where the last relation follows from the definition of the QCD scale parameter A: 



as = as{fi)= [^Sln^j . (19) 

Recalling that z = fir, Eq. (|l^) tells us that the series begins to diverge when the integral 
in Eq. ([T0|) is dominated by distances of order A~^. This feature is reminiscent of infrared 
renormalons. 

Within the context of perturbation theory, the best one can do is to truncate the series 
at the minimum term. We take the size of this term as an estimate of the inherent ambiguity 
in the perturbative calculation of the energy. Inserting Eq. (|I8]) into Eq. (|1^) yields 

/ 9 \ 3/2 



The analogous calculation for the state with quantum numbers n,l,m with wavefunctionQ 
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yields the result 
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(with a given by Eq. @). This is the main result of the paper. 

In a full QCD calculation, the ambiguity in the perturbative calculation of the energy, 
Eq. (p3D , would be cancelled by a related ambiguity in the nonperturbative contribution 
to the energy, resulting in an unambiguous expression. This cancellation between pertur- 
bative and nonperturbative ambiguities is in the same spirit as the analogous cancellation 
which occurs for infrared renormalons Thus we learn that quarkonia energies have a 
nonperturbative contribution which scales with naA like Eq. (0). However, it is impos- 
sible to calculate the coefficient of this nonperturbative contribution within the context of 
perturbation theory. 

The ambiguity in the perturbative calculation of the energy, Eq. (^), is nearly propor- 
tional to the probability density that the quarks in quarkonia are separated by a distance 
of A~\ which is given by (p)p ~ p^'+^[L^'_j'/_^(/))]^e~'', with p = 2/naA. The average 

separation for the quarks is roughly n'^a. The ambiguity increases as p decreases, i.e., as the 
size of the state increases, until p ~ 0{n), at which point it apparently begins to decrease. 
However, this value of p corresponds to n^a ~ A~^, where our perturbative approach is 
unreliable, and Eq. ( p3D should not be trusted. 

The ambiguity in the perturbative calculation of the energy, Eq. (^31) , is proportional 
to e"^/""^ rather than a power of aA (as would arise from an infrared renormalon) because 
the series for the energy correction, Eq. (pISD, diverges less rapidly than factorially. An 
approximate solution to Eq. (|T3|) is zq ^ kc/ In kc; thus the factor In'^ zo e~^°^'^ in Eq. (|T5|) is 
approximately \n^ke~'', which grows less rapidly than factorially, k\ ~ k^ e'^ . A similarly 
divergent series, resulting in a perturbative ambiguity proportional to e"'^'^^"'^)/^, appears in 
the calculation of soft-gluon resummation in the production of an object of mass Q from a 
hadron colhsion of total energy (r = Q'^/S) 

A nonperturbative contribution to quarkonia energies from the gluon condensate has 
been previously studied 0, ^. This contribution is associated with an infrared renormalon 
arising from emission and absorption of soft gluons, as shown in Fig. 2. The perturbative 
ambiguity from this renormalon is cancelled by a related ambiguity in the gluon condensate 
0. The gluon-condensate contribution to the quarkonia energies is 0, |^ 

AEni = ^rriQn^a^ < a.G^ > e^i (24) 



^We use the notation and conventions of Ref. The L^l^l^i are associated Laguerre polynomials. 
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Figure 2: A higher-order contribution to the quarkonium energy from the emission and 
absorption of gluons. The insertion of vacuum-polarization subdiagrams in the gluon prop- 
agator leads to an infrared renormalon associated with the gluon condensate. 



where e^/ are numbers of order unity. Since < agC^ > is of order A^, this contribution is 
proportional to (aA)^. Both the physical origin of Eq. and its functional dependence on 
aA show that it is a completely separate nonperturbative effect from the one considered in 
this paper. 

The gluon condensate has been considered as the leading nonperturbative contribution 
to quarkonia in various analyses [§, |l^, |lT[ |l3l. For states which are small compared 
with A~^, the gluon condensate, proportional to (aA)^, is formally more important than the 
nonperturbative effect discussed in this paper, proportional to e~^/""'^. However, as n^a ap- 
proaches A~^, the relative size of the gluon-condensate contribution and the nonperturbative 
contribution considered in this paper are unknown, since we do not know the normalization 
of the latter. It is therefore important to establish the size of this nonperturbative effect. 

In this paper we have shown that the calculation of quarkonia energies is divergent at 
large orders in perturbation theory. The associated ambiguity in the perturbative result, 
Eq. (p3D, is proportional to e"^/""^, where a is the Bohr radius, n is the principal quantum 
number, and A is the QCD scale parameter. This ambiguity is associated with a nonper- 
turbative contribution to the energy, coming from distances of order A~^ and greater, which 



scales with naA like Eq. (23). This is a separate nonperturbative effect from the gluon- 
condensate contribution to quarkonia energies. It is important to establish the size of this 
nonperturbative effect to judge its impact on quarkonia energies. 
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